Genome-Wide Mining of CULLIN E3 Ubiquitin Ligase Genes from Uncaria rhynchophylla.

CULLIN (CUL) protein is a subtype of E3 ubiquitin ligase that is involved in a variety of biological processes and responses to stress in plants. In Uncaria rhynchophylla, the CUL gene family has not been identified and its role in plant development, stress response and secondary metabolite synthesis has not been studied. In this study, 12 UrCUL gene members all contained the typical N-terminal domain and C-terminal domain identified from the U. rhynchophylla genome and were classified into four subfamilies based on the phylogenetic relationship with CULs in Arabidopsis thaliana. They were unevenly distributed on eight chromosomes but had a similar structural composition in the same subfamily, indicating that they were relatively conserved and potentially had similar gene functions. An interspecific and intraspecific collinearity analysis showed that fragment duplication played an important role in the evolution of the CUL gene family. The analysis of the cis-acting elements suggests that the UrCULs may play an important role in various biological processes, including the abscisic acid (ABA) response. To investigate this hypothesis, we treated the roots of U. rhynchophylla tissue-cultured seedlings with ABA. The expression pattern analysis showed that all the UrCUL genes were widely expressed in roots with various expression patterns. The co-expression association analysis of the UrCULs and key enzyme genes in the terpenoid indole alkaloid (TIA) synthesis pathway revealed the complex expression patterns of 12 UrCUL genes and some key TIA enzyme genes, especially UrCUL1, UrCUL1-likeA, UrCUL2-likeA and UrCUL2-likeB, which might be involved in the biosynthesis of TIAs. The results showed that the UrCULs were involved in the response to ABA hormones, providing important information for elucidating the function of UrCULs in U. rhynchophylla. The mining of UrCULs in the whole genome of U. rhynchophylla provided new information for understanding the CUL gene and its function in plant secondary metabolites, growth and development.

Identification and characterization of two O-methyltransferases involved in methylated 2-(2-phenethyl) chromones biosynthesis in agarwood.

The 2-(2-phenethyl)chromones (PECs) are the signature constituents responsible for the fragrance and pharmacological properties of agarwood. O-methyltransferases (OMTs) are necessary for the methylated PECs biosynthesis, however, there is little information known about OMTs in Aquilaria sinensis. In this study, we identified 29 OMT genes from the A. sinensis genome. Expression analysis showed they were differentially expressed in different tissues and responded to drill wounding. Comprehensive analysis of the gene expression and methylated PECs content revealed that the AsOMT2, AsOMT8, AsOMT11, AsOMT16, and AsOMT28 could potentially be involved in methylated PECs biosynthesis. The vitro enzyme assays and functional analysis in Nicotiana benthamiana demonstrated that AsOMT11 and AsOMT16 could methylate 6-hydroxy-2-(2-phenylethyl)chromone to form 6-methoxy-2-(2-phenylethyl) chromone. The transient overexpression experiment in Qi-Nan revealed that AsOMT11 and AsOMT16 could significantly promote the accumulation of three major methylated PECs. Our results provide candidate genes for the mass production of methylated PECs using synthetic biology.

Genome-wide identification of GATA transcription factor family and the effect of different light quality on the accumulation of terpenoid indole alkaloids in Uncaria rhynchophylla.

Uncaria rhynchophylla is an evergreen vine plant, belonging to the Rubiaceae family, that is rich in terpenoid indole alkaloids (TIAs) that have therapeutic effects on hypertension and Alzheimer's disease. GATA transcription factors (TF) are a class of transcription regulators that participate in the light response regulation, chlorophyll synthesis, and metabolism, with the capability to bind to GATA cis-acting elements in the promoter region of target genes. Currently the charactertics of GATA TFs in U. rhynchophylla and how different light qualities affect the expression of GATA and key enzyme genes, thereby affecting the changes in U. rhynchophylla alkaloids have not been investigated. In this study, 25 UrGATA genes belonging to four subgroups were identified based on genome-wide analysis. Intraspecific collinearity analysis revealed that only segmental duplications were identified among the UrGATA gene family. Collinearity analysis of GATA genes between U. rhynchophylla and four representative plant species, Arabidopsis thaliana, Oryza sativa, Coffea Canephora, and Catharanthus roseus was also performed. U. rhynchophylla seedlings grown in either red lights or under reduced light intensity had altered TIAs content after 21 days. Gene expression analysis reveal a complex pattern of expression from the 25 UrGATA genes as well as a number of key TIA enzyme genes. UrGATA7 and UrGATA8 were found to have similar expression profiles to key enzyme TIA genes in response to altered light treatments, implying that they may be involved in the regulation TIA content. In this research, we comprehensively analyzed the UrGATA TFs, and offered insight into the involvement of UrGATA TFs from U. rhynchophylla in TIAs biosynthesis.

Reference Genes Screening and Gene Expression Patterns Analysis Involved in Gelsenicine Biosynthesis under Different Hormone Treatments in Gelsemium elegans.

Reverse transcription quantitative polymerase chain reaction (RT-qPCR) is an accurate method for quantifying gene expression levels. Choosing appropriate reference genes to normalize the data is essential for reducing errors. Gelsemium elegans is a highly poisonous but important medicinal plant used for analgesic and anti-swelling purposes. Gelsenicine is one of the vital active ingredients, and its biosynthesis pathway remains to be determined. In this study, G. elegans leaf tissue with and without the application of one of four hormones (SA, MeJA, ETH, and ABA) known to affect gelsenicine synthesis, was analyzed using ten candidate reference genes. The gene stability was evaluated using GeNorm, NormFinder, BestKeeper, ∆CT, and RefFinder. The results showed that the optimal stable reference genes varied among the different treatments and that at least two reference genes were required for accurate quantification. The expression patterns of 15 genes related to the gelsenicine upstream biosynthesis pathway was determined by RT-qPCR using the relevant reference genes identified. Three genes 8-HGO, LAMT, and STR, were found to have a strong correlation with the amount of gelsenicine measured in the different samples. This research is the first study to examine the reference genes of G. elegans under different hormone treatments and will be useful for future molecular analyses of this medically important plant species.

Evaluation of Reference Genes for Normalizing RT-qPCR and Analysis of the Expression Patterns of WRKY1 Transcription Factor and Rhynchophylline Biosynthesis-Related Genes in Uncaria rhynchophylla.

Uncaria rhynchophylla (Miq.) Miq. ex Havil, a traditional medicinal herb, is enriched with several pharmacologically active terpenoid indole alkaloids (TIAs). At present, no method has been reported that can comprehensively select and evaluate the appropriate reference genes for gene expression analysis, especially the transcription factors and key enzyme genes involved in the biosynthesis pathway of TIAs in U. rhynchophylla. Reverse transcription quantitative PCR (RT-qPCR) is currently the most common method for detecting gene expression levels due to its high sensitivity, specificity, reproducibility, and ease of use. However, this methodology is dependent on selecting an optimal reference gene to accurately normalize the RT-qPCR results. Ten candidate reference genes, which are homologues of genes used in other plant species and are common reference genes, were used to evaluate the expression stability under three stress-related experimental treatments (methyl jasmonate, ethylene, and low temperature) using multiple stability analysis methodologies. The results showed that, among the candidate reference genes, S-adenosylmethionine decarboxylase (SAM) exhibited a higher expression stability under the experimental conditions tested. Using SAM as a reference gene, the expression profiles of 14 genes for key TIA enzymes and a WRKY1 transcription factor were examined under three experimental stress treatments that affect the accumulation of TIAs in U. rhynchophylla. The expression pattern of WRKY1 was similar to that of tryptophan decarboxylase (TDC) under ETH treatment. This research is the first to report the stability of reference genes in U. rhynchophylla and provides an important foundation for future gene expression analyses in U. rhynchophylla. The RT-qPCR results indicate that the expression of WRKY1 is similar to that of TDC under ETH treatment. It may coordinate the expression of TDC, providing a possible method to enhance alkaloid production in the future through synthetic biology.

Optimization of Isolation and Transformation of Protoplasts from Uncaria rhynchophylla and Its Application to Transient Gene Expression Analysis.

Protoplast-based engineering has become an important tool for basic plant molecular biology research and developing genome-edited crops. Uncaria rhynchophylla is a traditional Chinese medicinal plant with a variety of pharmaceutically important indole alkaloids. In this study, an optimized protocol for U. rhynchophylla protoplast isolation, purification, and transient gene expression was developed. The best protoplast separation protocol was found to be 0.8 M D-mannitol, 1.25% Cellulase R-10, and 0.6% Macerozyme R-10 enzymolysis for 5 h at 26 °C in the dark with constant oscillation at 40 rpm/min. The protoplast yield was as high as 1.5 × 107 protoplasts/g fresh weight, and the survival rate of protoplasts was greater than 90%. Furthermore, polyethylene glycol (PEG)-mediated transient transformation of U. rhynchophylla protoplasts was investigated by optimizing different crucial factors affecting transfection efficiency, including plasmid DNA amount, PEG concentration, and transfection duration. The U. rhynchophylla protoplast transfection rate was highest (71%) when protoplasts were transfected overnight at 24 °C with the 40 µg of plasmid DNA for 40 min in a solution containing 40% PEG. This highly efficient protoplast-based transient expression system was used for subcellular localization of transcription factor UrWRKY37. Finally, a dual-luciferase assay was used to detect a transcription factor promoter interaction by co-expressing UrWRKY37 with a UrTDC-promoter reporter plasmid. Taken together, our optimized protocols provide a foundation for future molecular studies of gene function and expression in U. rhynchophylla.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Siraitia siamensis.

WRKY transcription factors, as the largest gene family in higher plants, play an important role in various biological processes including growth and development, regulation of secondary metabolites, and stress response. In this study, we performed genome-wide identification and analysis of WRKY transcription factors in S. siamensis. A total of 59 SsWRKY genes were identified that were distributed on all 14 chromosomes, and these were classified into three major groups based on phylogenetic relationships. Each of these groups had similar conserved motifs and gene structures. We compared all the S. siamensis SsWRKY genes with WRKY genes identified from three diverse plant species, and the results implied that segmental duplication and tandem duplication play an important roles in the evolution processes of the WRKY gene family. Promoter region analysis revealed that SsWRKY genes included many cis-acting elements related to plant growth and development, phytohormone response, and both abiotic and biotic stress. Expression profiles originating from the transcriptome database showed expression patterns of these SsWRKY genes in four different tissues and revealed that most genes are expressed in plant roots. Fifteen SsWRKY genes with low-temperature response motifs were surveyed for their gene expression under cold stress, showing that most genes displayed continuous up-regulation during cold treatment. Our study provides a foundation for further study on the function and regulatory mechanism of the SsWRKY gene family.

Genome-Wide Identification of m6A Writers, Erasers and Readers in Poplar 84K.

N6-methyladenosine (m6A) RNA modification is a conserved mechanism to regulate gene expression that plays vital roles in the development of plants. However, the m6A RNA modification in forest trees remains limited. Here, we performed a complete analysis of m6A writers, erasers and readers in Poplar 84K, including gene location, gene structures, conserved motifs, phylogenetic relationships, promoter analysis, expression profiles and the homology modeling. We have identified 61 m6A pathway genes in Poplar 84K (Populus alba × Populus glandulosa), including 14 m6A writers, 14 m6A erasers and 33 m6A readers. Phylogenetic analysis indicated that the m6A writers and erasers were clustered into four groups and m6A readers were clustered into two groups. Promoter analysis showed that m6A pathway genes were mainly responsive to low oxygen followed by ABA and ethylene. The expression of the identified m6A pathway genes showed tissue-specific expression patterns in leaves, xylem, phloem and roots. Moreover, 17 genes were significantly up-regulated and 13 genes were significantly down-regulated in poplar overexpressing the transcription factor LBD15. Homology modeling and molecular docking results suggested that PagFIP37b was most likely to be regulated by LBD15, and the qPCRshowed that PagFIP37b were up-regulated in the LBD15-oe plants. The results provide insights that aid in the future elucidation of the functions of these m6A pathway genes and the epigenetic regulation mechanism of these genes in Poplar 84K.

The Research Progress of Taxol in Taxus.

BACKGROUND: Taxus is a valuable woody species with important medicinal value. The bark of Taxus can produce taxol, a natural antineoplastic drug that is widely used in the treatment of breast, ovarian and lung cancers. However, the low content of taxol in the bark of Taxus can not meet the growing clinical demands, so the current research aims at finding ways to increase taxol production. OBJECTIVE: In this review, the research progress of taxol including the factors affecting the taxol content, biosynthesis pathway of taxol, production of taxol in vitro and the application of multi-omics approaches in Taxus as well as future research prospects will be discussed. RESULTS: The taxol content is not only dependent on the species, age and tissues but is also affected by light, moisture levels, temperature, soil fertility and microbes. Most of the enzymes in the taxol biosynthesis pathway have been identified and characterized. Total chemical synthesis, semi-synthesis, plant cell culture and biosynthesis in endophytic fungi have been explored to product taxol. Multi-omics have been used to study Taxus and taxol. CONCLUSION: Further efforts in the identification of unknown enzymes in the taxol biosynthesis pathway, establishment of the genetic transformation system in Taxus and the regulatory mechanism of taxol biosynthesis and Taxus cell growth will play a significant role in improving the yield of taxol in Taxus cells and plants.

Selection and validation of reference genes for quantitative gene expression normalization in Taxus spp.

Quantitative real-time PCR (qRT-PCR) is commonly used to measure gene expression to further explore gene function, while suitable reference genes must be stably expressed under different experimental conditions to obtain accurate and reproducible data for relative quantification. Taxol or paclitaxel is an important anticancer compound mainly identified in Taxus spp. The molecular mechanism of the regulation of taxol biosynthesis is current research goal. However, in the case of Taxus spp., few reports were published on screening suitable reference genes as internal controls for qRT-PCR. Here, eight reference genes were selected as candidate reference genes for further study. Common statistical algorithms geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder were used to analyze the data from samples collected from a cell line of Taxus × media under various experimental conditions and from tissues of Taxus chinensis var. mairei. The expression patterns of TcMYC under salicylic acid treatment differed significantly, with the best and worst reference genes in the cell line. This study screened out suitable reference genes (GAPDH1 and SAND) under different treatments and tissues for the accurate and reliable normalization of the qRT-PCR expression data of Taxus spp. At the same time, this study will aid future research on taxol biosynthesis-related genes expression in Taxus spp., and can also be directly used to other related species.

The R2R3-MYB transcription factor family in Taxus chinensis: identification, characterization, expression profiling and posttranscriptional regulation analysis.

The MYB transcription factor family is one of the largest gene families playing regulatory roles in plant growth and development. The MYB family has been studied in a variety of plant species but has not been reported in Taxus chinensis. Here we identified 72 putative R2R3-MYB genes in T. chinensis using a comprehensive analysis. Sequence features, conversed domains and motifs were characterized. The phylogenetic analysis showed TcMYBs and AtMYBs were clustered into 36 subgroups, of which 24 subgroups included members from T. chinensis and Arabidopsis thaliana, while 12 subgroups were specific to one species. This suggests the conservation and specificity in structure and function of plant R2R3-MYBs. The expression of TcMYBs in various tissues and different ages of xylem were investigated. Additionally, miRNA-mediated posttranscriptional regulation analysis revealed that TcMYBs were the targets of miR858, miR159 and miR828, suggesting the posttranscriptional regulation of MYBs is highly conserved in plants. The results provide a basis for further study the role of TcMYBs in the regulation of secondary metabolites of T. chinensis.

Identification of the Genes Involved in Anthocyanin Biosynthesis and Accumulation in Taxus chinensis.

Taxus chinensis is a precious woody species with significant economic value. Anthocyanin as flavonoid derivatives plays a crucial role in plant biology and human health. However, the genes involved in anthocyanin biosynthesis have not been identified in T. chinensis. In this study, twenty-five genes involved in anthocyanin biosynthesis were identified, including chalcone synthase, chalcone isomerase, flavanone 3-hydroxylase, anthocyanidin synthase, flavonoid 3'-hydroxylase, flavonoid 3',5'-hydroxylase, dihydroflavonol 4-reductase, anthocyanidin reductase, and leucoanthocyanidin reductase. The conserved domains and phylogenetic relationships of these genes were characterized. The expression levels of these genes in different tissues and different ages of xylem were investigated. Additionally, the anthocyanin accumulation in xylem of different ages of T. chinensis was measured. The results showed the anthocyanin accumulation was correlated with the expression levels of dihydroflavonol 4-reductase, anthocyanidin synthase, flavonoid 3'-hydroxylase, and flavonoid 3',5'-hydroxylase. Our results provide a basis for studying the regulation of the biosynthetic pathway for anthocyanins and wood color formation in T. chinensis.

A comparative metabolomics analysis of the components of heartwood and sapwood in Taxus chinensis (Pilger) Rehd.

Taxus chinensis is a well-known gymnosperm with great ornamental and medicinal value. Its purple red brown heartwood (HW) has many attributes such as straight texture, high density, mechanical strength, rich elasticity and corrosion resistance that is highly prized commercially. T. chinensis sapwood (SW), in comparison, lacks these important traits. At present, little is known about the differences of metabolites between the SW and HW in T. chinensis. Widely targeted metabolic profiling was performed to analyze the metabolic profiles of HW and SW in T. chinensis using Liquid Chromatography-Electrospray Ionization-Mass Spectrometry (LC-EI-MS). A total of 607 metabolites were detected in HW and SW. Among them, 146 metabolites were significantly higher, and 167 metabolites significantly lower, in HW as compared to SW. These differential metabolites were mainly involved in metabolic pathways and biosynthesis of secondary metabolites, such as flavonoids, flavone and flavonol, phenylpropanoids and antibiotics. Moreover, 71 flavonoids and isoflavones were found to be significantly different between HW and SW. Our results show the difference of components between the HW and SW, which has potential significance to further elucidate the mechanism of HW color formation. The results will provide insight into the metabolites associated with wood color formation and useful information for understanding the metabolites associated with wood quality.

The genome of Populus alba x Populus tremula var. glandulosa clone 84K.

Poplar 84K (Populus alba x P. tremula var. glandulosa) is a fast-growing poplar hybrid. Originated in South Korea, this hybrid has been extensively cultivated in northern China. Due to the economic and ecological importance of this hybrid and high transformability, we now report the de novo sequencing and assembly of a male individual of poplar 84K using PacBio and Hi-C technologies. The final reference nuclear genome (747.5 Mb) has a contig N50 size of 1.99 Mb and a scaffold N50 size of 19.6 Mb. Complete chloroplast and mitochondrial genomes were also assembled from the sequencing data. Based on similarities to the genomes of P. alba var. pyramidalis and P. tremula, we were able to identify two subgenomes, representing 356 Mb from P. alba (subgenome A) and 354 Mb from P. tremula var. glandulosa (subgenome G). The phased assembly allowed us to detect the transcriptional bias between the two subgenomes, and we found that the subgenome from P. tremula displayed dominant expression in both 84K and another widely used hybrid, P. tremula x P. alba. This high-quality poplar 84K genome will be a valuable resource for poplar breeding and for molecular biology studies.

Transcriptomic Analysis of Betula halophila in Response to Salt Stress.

Soil salinization is a matter of concern worldwide. It can eventually lead to the desertification of land and severely damage local agricultural production and the ecological environment. Betula halophila is a tree with high salt tolerance, so it is of importance to understand and discover the salt responsive genes of B. halophila for breeding salinity resistant varieties of trees. However, there is no report on the transcriptome in response to salt stress in B. halophila. Using Illumina sequencing platform, approximately 460 M raw reads were generated and assembled into 117,091 unigenes. Among these unigenes, 64,551 unigenes (55.12%) were annotated with gene descriptions, while the other 44.88% were unknown. 168 up-regulated genes and 351 down-regulated genes were identified, respectively. These Differentially Expressed Genes (DEGs) involved in multiple pathways including the Salt Overly Sensitive (SOS) pathway, ion transport and uptake, antioxidant enzyme, ABA signal pathway and so on. The gene ontology (GO) enrichments suggested that the DEGs were mainly involved in a plant-type cell wall organization biological process, cell wall cellular component, and structural constituent of cell wall molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment showed that the top-four enriched pathways were 'Fatty acid elongation', 'Ribosome', 'Sphingolipid metabolism' and 'Flavonoid biosynthesis'. The expression patterns of sixteen DEGs were analyzed by qRT-PCR to verify the RNA-seq data. Among them, the transcription factor AT-Hook Motif Nuclear Localized gene and dehydrins might play an important role in response to salt stress in B. halophila. Our results provide an important gene resource to breed salt tolerant plants and useful information for further elucidation of the molecular mechanism of salt tolerance in B. halophila.

Genomewide analysis of the lateral organ boundaries domain gene family in Eucalyptus grandis reveals members that differentially impact secondary growth.

Lateral Organ Boundaries Domain (LBD) proteins are plant-specific transcription factors playing crucial roles in growth and development. However, the function of LBD proteins in Eucalyptus grandis remains largely unexplored. In this study, LBD genes in E. grandis were identified and characterized using bioinformatics approaches. Gene expression patterns in various tissues and the transcriptional responses of EgLBDs to exogenous hormones were determined by qRT-PCR. Functions of the selected EgLBDs were studied by ectopically overexpressing in a hybrid poplar (Populus alba × Populus glandulosa). Expression levels of genes in the transgenic plants were investigated by RNA-seq. Our results showed that there were forty-six EgLBD members in the E. grandis genome and three EgLBDs displayed xylem- (EgLBD29) or phloem-preferential expression (EgLBD22 and EgLBD37). Confocal microscopy indicated that EgLBD22, EgLBD29 and EgLBD37 were localized to the nucleus. Furthermore, we found that EgLBD22, EgLBD29 and EgLBD37 were responsive to the treatments of indol-3-acetic acid and gibberellic acid. More importantly, we demonstrated EgLBDs exerted different influences on secondary growth. Namely, 35S::EgLBD37 led to significantly increased secondary xylem, 35S::EgLBD29 led to greatly increased phloem fibre production, and 35S::EgLBD22 showed no obvious effects. We revealed that key genes related to gibberellin, ethylene and auxin signalling pathway as well as cell expansion were significantly up- or down-regulated in transgenic plants. Our new findings suggest that LBD genes in E. grandis play important roles in secondary growth. This provides new mechanisms to increase wood or fibre production.

Genome-wide identification and characterization of the SPL gene family in Ziziphus jujuba.

SQUAMOSA Promoter-Binding Protein-Likes (SPLs) are plant specific transcription factors playing important roles in plant growth and development. The SPL gene family has been studied in various plant species; however, there is no report about SPLs in Zizyphus jujuba. In this study, we identified 18 putative ZjSPL genes in Z. jujuba using a genome-wide analysis. Sequence features, gene structures, conserved domains and motifs were analyzed. The phylogenetic relationships of SPLs in Z. jujuba and A. thaliana were revealed. A total of 5 pairs of ZjSPLs were identified, suggesting the importance of gene duplication in SPL gene expansion in Z. jujuba. In addition, 11 of the 18 ZjSPLs, belonging to G1, G2 and G5 subgroups, were found to be targets of miR156, suggesting the conservation of miR156-mediated posttranscriptional regulation in plants. Expression analysis revealed that eight ZjSPL genes were responsive to the infection of witches'-broom phytoplasma. Our results provide a basis for the further elucidation of the biological function of ZjSPLs and their regulation in witches'-broom disease.

Genome-Wide Identification and Analysis of MicroRNAs Involved in Witches'-Broom Phytoplasma Response in Ziziphus jujuba.

MicroRNAs (miRNAs) play an important role in responding to biotic and abiotic stresses in plants. Jujube witches'-broom a phytoplasma disease of Ziziphus jujuba is prevalent in China and is a serious problem to the industry. However, the molecular mechanism of the disease is poorly understood. In this study, genome-wide identification and analysis of microRNAs in response to witches'-broom was performed. A total of 85 conserved miRNA unique sequences belonging to 32 miRNA families and 24 novel miRNA unique sequences, including their complementary miRNA* strands were identified from small RNA libraries derived from a uninfected and witches'-broom infected Z. jujuba plant. Differentially expressed miRNAs associated with Jujube witches'-broom disease were investigated between the two libraries, and 12 up-regulated miRNAs and 10 down- regulated miRNAs identified with more than 2 fold changes. Additionally, 40 target genes of 85 conserved miRNAs and 49 target genes of 24 novel miRNAs were predicted and their putative functions assigned. Using the modified 5'-RACE method, we confirmed that SPL and MYB were cleaved by miR156 and miR159, respectively. Our results provide insight into the molecular mechanisms of witches'-broom disease in Z. jujuba.

Identification of novel and conserved microRNAs in Panax notoginseng roots by high-throughput sequencing.

BACKGROUND: MicroRNAs (miRNAs) are small, non-coding RNAs that are important regulators of gene expression, and play major roles in plant development and their response to the environment. Root extracts from Panax notoginseng contain triterpene saponins as their principal bioactive constituent, and demonstrate medicinal properties. To investigate the novel and conserved miRNAs in P. notoginseng, three small RNA libraries constructed from 1-, 2-, and 3-year-old roots in which root saponin levels vary underwent high-throughput sequencing. METHODS: P. notoginseng roots, purified from 1-, 2-, and 3-year-old roots, were extracted for RNA, respectively. Three small libraries were constructed and subjected to next generation sequencing. RESULTS: Sequencing of the three libraries generated 67,217,124 clean reads from P. notoginseng roots. A total of 316 conserved miRNAs (belonging to 67 miRNA families and one unclassified family) and 52 novel miRNAs were identified. MIR156 and MIR166 were the largest miRNA families, while miR156i and miR156g showed the highest abundance of miRNA species. Potential miRNA target genes were predicted and annotated using Cluster of Orthologous Groups, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes. Comparing these miRNAs between root samples revealed 33 that were differentially expressed between 2- and 1-year-old roots (8 increased, 25 decreased), 27 differentially expressed between 3- and 1-year-old roots (7 increased, 20 decreased), and 29 differentially expressed between 3- and 2-year-old roots (8 increased, 21 decreased). Two significantly differentially expressed miRNAs and four miRNAs predicted to target genes involved in the terpenoid backbone biosynthesis pathway were selected and validated by quantitative reverse transcription PCR. Furthermore, the expression patterns of these six miRNAs were analyzed in P. notoginseng roots, stems, and leaves at different developmental stages. CONCLUSIONS: This study identified a large number of P. notoginseng miRNAs and their target genes, functional annotations, and gene expression patterns. It provides the first known miRNA profiles of the P. notoginseng root development cycle.

Comparative analysis of the Dicer-like gene family reveals loss of miR162 target site in SmDCL1 from Salvia miltiorrhiza.

DCL1, the core component for miRNA biogenesis, is itself regulated by miR162 in Arabidopsis. MiRNA-mediated feedback regulation of AtDCL1 is important to maintain the proper level of DCL1 transcripts. However, it is unknown whether the miRNA-mediated regulation of DCL1 is conserved among plants. We analyzed the SmDCL gene family in Salvia miltiorrhiza, an emerging model plant for Traditional Chinese Medicine (TCM) studies, using a comprehensive approach integrating genome-wide prediction, molecular cloning, gene expression profiling, and posttranscriptional regulation analysis. A total of five SmDCLs were identified. Comparative analysis of SmDCLs and AtDCLs showed an apparent enlargement of SmDCL introns in S. miltiorrhiza. The absence of miR162 in S. miltiorrhiza and the loss of miR162 target site in SmDCL1 were unexpectedly found. Further analysis showed that the miR162 target site was not present in DCL1 from ancient plants and was gained during plant evolution. The gained miR162 target site might be lost in a few modern plants through nucleotide mutations. Our results provide evidence for the gain and loss of miR162 and its target sites in Dicer-like genes during evolution. The data is useful for understanding the evolution of miRNA-mediated feedback regulation of DCLs in plants.